function [RFE,convKernels,T,v]=timeInvariantSingleCell(x,y,numLags,pen)
if nargin<4
    pen=0;
end
[N,bins]=size(x);
pEE=zeros(bins,1);
for b=1:bins
    %     ndx=[1:b-1,b+1:bins];
    [convKernels(:,b),T(:,:,b),v(:,b),RFE(b)]=getBoundedWiener(x(:,b),y(:,b),numLags,pen);
    %     pEE(b)=getErr(aStationary,x(:,b),y(:,b));
end